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In this paper we calculate the scale anomaly for a quantum field theoretic 2D non-relativistic Bose 
gas with contact interactions using Fujikawa’s method, both in vacuum and in many-body systems. 
The use of path integrals for these problems is novel and motivated by a recently developed path- 
integral framework for addressing questions about scaling in these systems. A natural class of 
regulators is found that produces the correct value of the anomaly traditionally calculated via other 
methods, e.g., diagrammatically via the P function. 
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I. INTRODUCTION 

The use of Fujikawa’s method in particle physics is well known and is now standard in textbooks [l|. It was originally 
developed to understand the chiral anomaly 0 but has since been extended to other cases, including the relativistic 
scale anomaly Q. However, as far as we are aware, it has not been used before for non-relativistic physics. There 
are currently reasons to embark in such calculations. Non-relativistic anomalies have been studied since the seminal 
paper by R. Jackiw Q, mostly using canonical methods, not Fujikawa’s^. Interest in these anomalies has intensified 
in the study of ultracold 2D g^es (THI^. with the work by J. Hofmann on anomalies of trapped 2D Fermion gases 
being of particular relevance [l6l| . Despite all this activity, there are still questions about anomalies and their impact 
in such systems that need to be answered El- A path-integral Fujikawa approach to study anomalies in systems 
with an SO{2, 1) classical symmetry, mainly in the context of 2D diluted gases, has been recently proposed in 
While this approach provides a nice picture of the structure of anomalies in many-body systems, the calculation of 
the Fujikawa Jacobian is crucial in order for this framework to also provide a practical scheme that will help us better 
understand the role of anomalies in lower-dimensional physics. We present here our first results of the Jacobian 
calculation for 2D complex fields with contact interactions in the case of constant background fields. 

Within the path-integral formulation, anomalies result from the presence of Jacobians due to the non-invariance 
of the measure under symmetry transformations. These Jacobians are functional determinants and need to be 
regularized. For the chiral anomaly, all regulators lead to a finite result, whereas for the relativistic scale anomaly 
an infinite piece remains that is present even if the same regulator is used in the free theory, so this piece can be 
subtracted if the free theory is taken to be non-anomalous [l^. The non-relativistic scale anomaly is similar to 
the relativistic case in this respect. However, unlike the latter, space and time are treated on unequal footing in 
the former. Indeed, traditionally, for both the relativistic chiral and scale anomalies, one goes into Euclidean space 
where the Lagrangian kinetic operator is Hermitian. In this Euclidean space one can work with functions of a single 
variable (the 4 momentum squared) that is positive semi-definite in all directions. In contrast, for the non-relativistic 
case the Lagrangian operator is Hermitian in real time (“Minkowski space”). Due to the asymmetry between space 
and time, one is stuck with lu and rather than a single k^, making the task considerably more difficult, which may 
be a reason for why this problem has not been addressed before using Eujikawa’s method. 

The structure of this paper is as follows: we give a brief introduction to Eujikawa’s method, after which we review the 
essential technical details for the system that will be considered here. We then proceed with the Jacobian calculation 
for zero and finite temperature. Conclusions and comments end the paper. 


*Electronic address: cllin@uh.edu Electronic address: cordonez@central.uh.edu 
^ See @,@1 and references therein. 
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II. FUJIKAWA’S DERIVATION 

The derivation of the anomaly via Fujikawa’s method presented here follows closely the path-integral derivation of 
the Ward identities, but now the Jacobian of the symmetry transformation is taken into account. Indeed, anomalies 
represent a breakdown of the Ward identities, and it is precisely the Jacobian that invalidates the identities. For 
simplicity we will demonstrate the derivation for a scalar field theory without sources: the generalization to other 
(multiple) fields is straightforward. With a change of variables given by (/''(x) = (j){x) + r]S(j){x): 
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Had we added a source term f d'^x J{x)4){x), the equation would read: 
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Differentiation w.r.t. to J{xi) n times and setting J = 0 would create contact terms: 
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Eq. (5) without the Jacobian contribution is the traditional Ward identity at zero temperature, in vacuum, presented 
in most textbooks [ 2 ^. In our case, we only need the Jacobian of the infinitesimal transformation by itself in order 
to compute the RHS of Eq. (3) and compare our results with the literature for both the zero-temperature and the 
finite-temperature case. For the latter, we will work within the framework of reference [Tsj l. for which a detailed 
calculation is mandatory. 
























3 


III. CONTACT INTERACTION 

The Schrodinger Lagrangian density for bosons with contact interaction in 2D is given by: 

£ = ^ 

which is the 2-body interaction with a V{x — y) = g 6^ {x — y) potential: 


( 6 ) 
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( 7 ) 


The action corresponding to this Lagrangian is scale-invariant. This can be readily seen by noting that in D = 2, the 
coupling g has no dimensions in units of length (with Ti = m = 1). Therefore Eq. (3) applies. 


IV. SCALE TRANSFORMATION 


Under a non-relativistic dilation transformation 


Setting A = 1 -f 77 for infinitesimal 77 : 


x' = Xx, 
t' = X^t, 


Sx = gx, 

St = 2gt, 

Sip = gOip{t, x) = gSip, 

Sip* = g6ip*(t,x) = gSip*, 

9 = ( -y - x-V- 2tdt ) . 


where D = d — 1 is the spatial dimension^. In this paper we will set D = 2. Therefore^: 
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Note that unlike translations, for dilations both conventions — Sip = g (^—1 — x ■'V — 2tdt^ ip(t,x) or Sip = 

77 ^1 -f a; • V -I- 2tdt^ 'P{t,x) — leading to currents of opposite sign, are widely used. We’ve adopted the former, 
which leads to a dilation charge of [1^ : 


D = y cfx X ■ j — 2tH, 
j = (ip'^Vip - Vt/^V) • 


( 11 ) 


^ We used S^p {S'lp*) for the infinitesimal change in -ip (V’*)? set Sip = rjSip to make the notation consistent with Eq. (2). 
^ Sometimes we write x = (xo,x) = (t,x) for notational convenience. 
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V. FUJIKAWA CALCULATION: SET UP 


The generalization of the scalar case to our Lagrangian is straightforward: 
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where we’ve used det A = 4 Comparison with Eq. (3) makes the generalization clear: 
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This expression is singular so needs to be regularized. This is done by expanding S^(x — y)l 2 , where I 2 = 
using the eigenbasis <j)n of a Hermitian operator M: 

d^(x-y)l2 = '^(l)„{xo,x)(l)i{yo,y). 

71 

Inserting a regulator that’s a function of M 
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with the property that R{0) = 1 so that at the end of the calculation we send A ^ 00 and lim R(^) = 1 
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{[M] = [A^]). The idea is to choose R such that large eigenvalues are suppressed giving a convergent sum: 

5\{x-x)l2 (t>n{.Xo,x)(l)\{xQ,X). (16) 

Once the Hermitian operator M has been selected, then the sum over n in Eq. (15) gives: 
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so that 


tr [5|j(a: - y)l 2 \ = tr 




(18) 


For the class of regulators defined by Eqs. (19) and (20) in the next section - a consistent choice for the untrapped 
system - it will be shown that the non-trivial contribution to Eq. (18) will have an even integrand in both w and k 
(Eq. (28)). This means that the derivative terms in 9 will give a null contribution and the only term that will survive 
is given by Eq. (18). To see this, one should take the space-time derivatives in Eq. (21), and then set x = y; these 
terms will give odd contributions to the integrand in (w, k) space when multiplied by the even terms of Eq. (28). 

Therefore, Eq. (18) is the expression we aim to calculate. 


4 Tr includes both functional and matrix indices; tr only refers to the 2x2 matrix indices in Eqs. (12) and (13). 
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VI. FUJIKAWA CALCULATION: MODE EXPANSION 
A. Zero Temperature 


For our Hermitian matrix we will choose 
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-9'^ 
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where the helds in M are constant background fields, and ^ (x^ 




is the quadratic Lagrangian resulting 


from a saddle point expansion of the action about the background field ip, and x is the shift of the original field from 
ip. We’ve included a chemical potential g for many-body physics that explicitly breaks scale-invariance, but we can 
always set g = 0 and as we will demonstrate, the inclusion of /r has no effect on the anomaly. 


For our regulating function R we will choose 
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which clearly satisifies i?(0) = 1. Plugging in this regulator into Eq. (17) and Fourier expanding 5^{x — y) gives: 
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We will now take {xo,x) = {yo,y)', will make a change of variables Cj = -^ and k = and then replace the tildes 
since they are dummy indices (oj and k are now dimensionless): 
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For notational convenience we will write the above expression as: 
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To evaluate the inverse in Eq. (23), we will use the identity for matrix inverses {D + B) ^ = D ^ — {D ^B)D ^ + 
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Note that the ie makes D± invertible. 
So Eq. (22) becomes: 


tr [5«(0)/2] = ±A^J^J ^ tr {Dg^ - {Dg^B)Dg^ + [Dg^ B){Dg^ B)Dg^) , (26) 

where we terminated the series at two powers of B, since each additional power of B produces a ^ that the A'^ 
prefactor can’t offset. 

The first term in the series, Dg"^, when doing the integral over w, is independent of the coupling: 
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Therefore this term is also contained in the free-case, which we take to be anomaly-free. So we subtract this term 
when calculating the anomaly. The next term {Dg^B)Dg^ has no diagonal elements, so is traceless. 


The only term to calculate is the B){D^^B)D^^ term which produces: 
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Now A± in Eq. (24) can be safely taken to ±1 (A —>■ oo). Eor both ± cases, the result is the same: 
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Plugging this into Eq. (13) gives: 
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This can be compared with [2^ (for the case of constant background fields) by making the replacement g — 


Because both R {^) = (l ± ^ work as regulators, any linear combination such that their coefficients add to one 

works. Eor example: 
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We have also verified that the following regulators work: 


R 




(33) 


B. Many-Body 


Under the formalism developed in [T^®: 


2£-2P = tT[6HiO)l2]- (34) 

However, here the anomalous term is evaluated in Euclidean space using the finite temperature rules. Eq. (28) with 
plus chosen is: 
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When going to finite temperature, the difference is that we have —dr instead of idt- The effect is to replace ui with iuj 
in Eq. (35). That is, had we started directly in Euclidean space, we would still get Eq. (35), but with u) replaced by 
iuj stemming from —dr replacing idt in our regulator. The second change is that the integral is a sum since the modes 
are discrete, with a /3 factor resulting from writing the delta function as 6 ^{x — y) = I ( 277 )^^ ^ 

n 

So a sum replaces the integral. So had we started directly in Euclidean space, we would have a sum over frequencies 
instead of an integral: 


tr [5fl(0)/2] 


27r /3 
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where ujn = for bosons. The ie no longer matters and the summation is standard: 


tr [(5 k(0)/2] = 



which in the limit of large A gives: 


tr [i5i?(0)/2] 
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^ The factor of I3A in the denominator of Eq. (52) in [T 3 | cancels the Euclidean version of the factor Jq dr f d^x from Eq. (12) in this 
paper, since for constant background fields our class of regulators gives a constant value for tr ((5jj(0)/2)- The time and spatial derivatives 
in Os and 0 in paper [T 3 l give no contributions in this case (untrapped) as explained here. Notice Eq. (34) is a 2 x 2 version of Eq. (52) 

in Ei. 
















So plugging this into Eq. (34) 


2£-2P=- 


47r 


(40) 


which agrees with pHI ] with g ^ 2g. For the finite temperature case, there is some ambiguity in the continuation to 
Euclidean space that affects the sign, where A+ leads to the correct sign, and A- leads to the negative sign. We take 
the view that the zero-temperature limit must reproduce the vacuum result. 


VII. CONCLUSION 

Fujikawa’s path-integral method has been applied to the Schrodinger Lagrangian to describe anomalies for 2D non- 
relativistic, SO{2, 1) scale-invariant complex bosons with contact interactions. A class of natural regulators was 
identified that gives results consistent with those in the literature, obtained with other methods, in both zero and 
finite-temperature cases [13, This work was motivated by the recent formulation of Fujikawa’s approach to analyze 
the anomaly structure for 2D gases with SO{2, 1) classical symmetry (and other systems with such symmetry) [l8|, 
which is relevant in the study of ultracold 2D trapped gases [l3]- It was important, therefore, that we made contact 
with established work using other techniques. Further work is needed for a deeper understanding of this method and 
its possible applications. In particular, heat kernel techniques will be used to investigate trapped systems. Work on 
these issues is in progress [2^. 
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